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Abstract 


The  two  stage  life  test  of  Bulgren  and  Hewett  (1973)  for  the 

v 

mean  lifetime,  0 , of  an  exponentially  distributed  lifetime  is 

^ N r H 

essentially  a two  stage  version  of  type  II  censoring.  Since  the  (If-  1 - 

' — '■  " _ ' , V' 

test  decision  only  determines  whether  0 is  above  or  below  0n  , 

k-  \ 

it  may  be  of  interest  to  estimate  0 following  the  decision  using 
the  test  data.  Epstein  (1960a)  found  confidence  intervals  for 
under  type  II  censoring.  This  report  shows  that  Epstein's  intervals 
may  be  modified  slightly  and  used  to  provide  confidence  intervals 
for  0 following  the  two  stage  test.  The  resulting  confidence 


k, 


intervals  are  shown  to  be  conservative* 
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Confidence  Intervals  for  an  Exponential 
Parameter  from  a Two  Stage  Life  Test 


I.  Introduction. 

In  the  theory  of  reliability  and  life  testing,  a 
lifetime  x is  often  assumed  to  have  an  exponential 
distribution.  Thus  the  probability  density  function  of 
t may  be  expressed  as 


f(x;e)  = 


l e-x/e 
W 


if 


> o 


(l.D 


elsewhere 


where  0 > 0 is  the  expected  lifetime. 

In  a life  test  with  type  II  censoring,  n items 
are  placed  on  test  and  observed  for  failures  until  the 

failure  occurs.  Properties  of  this  test  were  given 
by  Epstein  and  Sobel  (1953) . Bulgren  and  Hewett  (1973) 
developed  a two  stage  version  of  a life  test  with  type  II 
censoring  to  test  Hq:  0 ^ eQ  vs.  H^:  0 < 0Q  . In  the 
test  n items  are  placed  on  test  and  two  integer  values, 
r^  and  are  chosen  such  that  r^  + ^ < n . In  the 

first  stage, testing  proceeds  (with  or  without  the  replacement 
of  failed  items)  until  r^  failures  occur.  At  this  time 
a decision  is  made  to  accept  or  reject  Hq  , or  to  continue 
on  to  the  second  stage  where  additional  failures  are 


1 


2 


observed  and  a final  decision  is  made.  The  decisions  are 

based  on  the  values  of  T and  T , the  total  times  on  test 

rl  r3 

at  the  times  of  the  and  failures  repectively, 

with  r^  = r^  + . The  test  is  depicted  in  figure  1.1. 

It  is  possible  that  additional  information  about 
the  true  value  of  0 may  be  needed  following  the  test 
decision.  Thus  it  may  be  of  interest  to  estimate  9 
using  the  data  from  the  life  test.  In  this  report  we 
derive  confidence  intervals  for  8 using  data  from  the 
two  stage  life  test. 


Stage 

1 


reject  Hq  j continue  to  stage  2j  accept  Hf 


2T 


Stage 

2 


reject  H. 


i accept  H. 
_i 9 


2T 


Figure  1.1.  The  two  stage  life  test. 


II.  A Two  Sided  Confidence  Interval  for  0 

Confidence  intervals  for  0 using  type  II  censoring 
were  given  by  Epstein  (1960a)  and  Epstein  (1960b).  A two 


3 


sided  interval  with  confidence  coefficient 


(1 


a)  is 

(2.1) 


where  is  the  total  time  on  test  at  the  rth  failure.  If 

tf , t2 t . . . , t represent  the  ordered  failure  times,  then 


r 

Z t.  + (n  - r)t  (without  replacement) 

i=l  1 r 

nt  (with  replacement) 

A two  sided  confidence  interval  for  type  II 
censoring  in  a single  stage  is  given  by  (2.1).  This 
suggests  the  possibility  that  the  single  stage  interval 
could  be  used  at  either  stage  of  the  two-stage  test.  More 
formally,  we  conjecture  that  the  following  rule  provides  a 
lOO(l-a')  percent  confidence  interval  for  0. 

When  the  decision  is  made  at  stage  1,  use 


4 


When  the  decision  is  made  at  stage  2,  use  (2.2) 


2T  2T 


We  shall  refer  to  the  first  interval  in  ( 2.2  ) as 

1^  and  the  second  interval  as  In  general  we  will  refer 

to  the  single  interval  generated  by  rule  ( 2.2  ) as  I.  To 
establish  that  I is  a confidence  interval  for  9 with  a 
coefficient  at  least  1-u',  we  would  have  to  show  that 
P(0el)  _>  1 " <*'•  Unfortunately,  it  will  be  shown  that 
P(0el)  < 1 - a'  for  some  0.  However,  we  will  also  show 
that,  for  a suitably  chosen  a' , the  confidence  interval 
given  by  ( 2.1  ) will  have  a confidence  coefficient  no  less 
than  1 - a. 

First,  define  event  D = (the  test  decision  is  made 
at  stage  1).  Then  5 will  represent  continuation  to  stage  2. 
We  may  then  write 

P(0el)  = PfDjOel^  + P(D,0eI2)  = p1  + P2>  say. 

2 

If  we  recall  that  2T  /9  ^ x->  • it  is  seen  that  P(0el1)  = 

rl  z 1 1 

1 - a'.  Also,  PfOel^  = PfD^el^  + P(D,0eI1>  = + 

P(D,0eI^).  Combining  these  results  gives 


4-H> 


. £ . t 
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P(6el)  = 1 - <x'  - (P (D#6e71)  - P(5f0el2)) 

= 1 - a'  - p*,  say. 

To  establish  the  validity  of  the  interval  in  ( 2.1  ) it 

would  be  necessary  that  p*  = P(D,0eI^)  - P(D,0eI2)  < 0. 

That  this  is  not  always  the  case  will  be  shown  later.  We 

first  introduce  some  notational  definitions  which  should 

make  subsequent  material  more  readable. 

Define  the  transformations  Y = 2T  /0  and 

rl 

Z = 2T  / 0 so  that  Y v x?  and  Z % • Under  these 

r 3 2r^  A2r^ 

transformations,  p*  is  equal  to 


6 


R 


3 


Thus,  under  the  transformations,  R^  , R^  and  R^  correspond  to 
1^,  the  continuation  region,  and  1^  respectively.  This 
permits  us  to  write  p*  simply  as  P(YcR2,YeR  ) - 
P(YcR2,ZrR3)  . 

Making  use  of  the  endpoints  of  intervals  R^ , R2 
and  R^ , we  can  make  some  general  observations  about 
P(0el).  We  consider  those  intervals  designated  by  numbers  in 
parentheses  in  Figure  2.1. 
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Figure  2.1.  A partition  of  the  parameter  space, 


In  interval  (1) : 


Vl 


0 < 


2 ,2t  3 
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In  this  interval. 


2 2 

nr1  > V > Xa ' 

— 2 r — 9 r 

2 ' 3 2 * ^ ^ 


which  implies  that  P(YeR2,ZeR3)  = P (YeR^ ,YeR^ ) = 0.  Thus, 
in  interval  (1)  we  have  P(8el)  = 1 - a' . 

In  interval  ( 2) : 


0nd  6 d 

< 0 < u 1 


xg 1 _ 

2 ' ^ r 3 


xj. 

I“'2rl 


Since 


Vi 


’ 4 2r 
2 '1 


in  this  interval,  we  have  P ( YeR^  rY'e  R^)  = 0.  This  gives  the 
result  that  P(0eZ)  _>  1 - a'  for  all  0 in  interval  (2). 

In  interval  ( 3)  : 


6 > 


60d2 

~T~ 

x a' 

' 2ri 


Here 


0Od2  „ 2 

q — X n 1 

H 9 r 

1 2 ' 1 
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implies  PCYtR^YeR^)  = 0.  So  again,  for  0 in  the  interval 
(3)  , P(0el)  > 1 - a'. 

In  interval  ( 4 ) : 


We  note  that,  if 


(d2/d1)  > 


SL.  2r 

2 ' 1 


/x,  «•  , 
1_T"'2  1 


the  sense  of  the  inequalities  will  be  reversed  and  interval 
(4)  will  be  empty.  When  it  is  non-empty  we  have  R^  0 R^  and 

P ( 6e  I ) = l-a'-P  (YeR2)  + P(YeR2,ZeR3) 

= l-a’-P(YeR2) (1-P(Z£R3!YeR2) ) 

= l-a'-P  (D)  (1-P(0el2  | D)  ) . 


Since  P(D)  > 0 and  P(0el2|5)  < 1,  it  follows  that 
P ( 0£ I ) < 1 - a'  for  0 in  interval  (4)  when  it  is  non-empty. 
This  shows  that  the  confidence  interval  proposed  by  (2.1) 
is  invalid  in  the  sense  that,  for  some  0,  P(0el)  < 1 - a' . 

For  those  regions  not  included  in  these  four 
intervals  a general  pattern  of  the  behavior  of  P(0el)  as  a 
function  of  0 was  observed  by  computing  P(0el)  on  the 
computer  for  various  test  parameters  given  by  Bulgren  and 


3 


Hewett  (1973).  The  general  pattern  which  was  observed  is 
depicted  in  Figure  2.2.  It  can  be  seen  from  Figure  2.2 
that  if  a'  were  chosen  sufficiently  small,  it  would  be 
possible  to  keep  P(9el)  _>  1 - a,  where  1 - a is  the  desired 
confidence  coefficient.  An  appropriate  value  of  a'  can  be 
found  if  we  find  an  upper  bound  on  p*  which  can  be  solved 
for  a'.  Recall  that  p*  is  just  F ( YeR2 , YeR^ ) - P ( Ye  R^ , Ze  R^)  . 
The  difficulty  in  finding  a bound  on  p*  is  finding  a closed 
expression  for  PfYeR^/ZeR^)  which  is  also  simple  enough  to 
permit  the  solution  for  a'  • Bulgren  and  Hewett  (1973)  give 
the  joint  density  of  Y and  Z as 


( z-y) 


r2“l  rL-l  -z/2 


ri  + r2 

2 <T(  rx)r(r2) 


0 <y  < z 0 < z<°° 


f(y,z)  = 


elsewhere . 


Then 


60d2  V 


2-'2r3 


P(YeR2,ZeR3)  = 


f (y , z) dzdy  , 


Vi 


X a, 





-i 


Vl 


eodi 


S0d2 


2 2 
(„  i i 

^ p V-  2r 
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eodi 

“2 
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Figure  2.2.  Behavior  of  P(0el)  when  I is  from 

(2.1) 

and  it  is  apparent  that  the  closed,  simple  expression  we 
seek  is  not  possible.  In  lieu  of  evaluating  this  joint 
probability  directly,  we  will  simplify  it  using  a Bonferrom 
bound.  A lower  bound  on  P (Yel^ , ZeR^)  , using  the  Bonferroni 
inequality,  is 

P(YeR2,ZeR3)  > P(YeR2)  + F(ZeR3)  - 1 
= P(YeR2)  + ( 1-a ' ) - 1 
= P(Y£R2)  - a1. 


Thus 


p*  < P(Y€R1,YcR2)  - P(YeR2)  + a'.  (2.3) 


11 


This  upper  bound  on  p*  may  be  further  simplified 
with  some  observations  about  the  relationship  of  R^  to  R2. 
With  respect  to  0,  interval  R1  is  fixed.  On  the  other  hand, 
^2  increases  in  both  size  and  location  as  0 decreases. 

When  0 = °°,  R2  is  the  single  point  {0}.  When  0 is  such 
that  0gd^/0  coincides  with 


/ 


it  could  be  that  R2S  f this  occurs,  the  upper  bound 

in  ( 2.3  ) achieves  its  maximum  value  a'.  Since  P(0el)  = 

1 - a'  - p*,  the  Bonferroni  inequality  gives  P(0eI)  > 

1 - 2a'.  Obviously,  if  we  want  P(0el)  _>  1 - a,  a'  is 
given  the  value  a/2.  This  establishes  the  following 
result. 


Theorem  2.1.  Let  r^  and  r2  be  the  number  of 

failures  required  for  each  stage  of  the  two-stage 

exponential  life  test,  where  r.  = r.  + r_.  If  T and  T 

3 i 2 rj  r3 

represent  the  total  time  on  test  at  the  r/h  and  r,th 

failures  respectively,  then  the  following  rule  gives  a 

two-sided  confidence  interval  for  the  exponential  mean  0, 

with  a confidence  coefficient  at  least  1 - a,  at  the  time 

of  decision  for  the  two-stage  test.  If  a decision  is 

made  at  the  first  stage,  use 


12 


If  a decision  is  made  at  the  second  stage,  use 


2T  2T 


The  wide  applicability  of  the  Bonferroni  inequality 
suggests  that  the  bound  in  ( 2.3  ) may  be  quite  crude, 
leading  to  overly  conservative  confidence  intervals.  An 
alternative  approach  will  be  employed  in  searching  for  a 
less  conservative  confidence  interval.  We  also  will  admit 
the  possibility  that  P(Qel)  < 1 - a,  for  some  0,  if  the 
difference  (l-ot)-P(6 el)  can  be  shown  to  be  small.  This 
would  result  in  approximate  confidence  intervals. 

Let 


Then 


A"  P(YeR2)  P (ZtR3)  - P(YeR2,  ZeR3)  . 


p*  = P(YeR1,YeR2)  - P ( Ye R2) P ( Ze R3)  + A 


= P(YeR1,YeR2)  - P(YeR2) (1-a')  + A.  (2.4) 


r— itwtwihw  limj,  „ 
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If  A = 0 (this  would  be  true  if  Y and  Z were  independent) 
p*  could  be  bounded  above  without  use  of  the  Bonferroni 
inequality.  The  penalty  for  assuming  A = 0,  when,  in  fact, 
A ^ 0 could  be  a confidence  interval  where  P(0el)  < 1 - a 
for  some  0.  The  difference  <!l-a)  - P(0el)  would  be  at  most 
A,  and  the  Bonferroni  inequality  can  provide  an  upper  bound 
on  A . 

P(YeR1,ZeR3)  > P(YeR2>  + P(ZeR3)  - 1 
= P(YeR2)  - a'  . 


Thus 


A < P(YeR2)  (1-a' ) - P(YeR2)  + a' 

= a' (l-P(YeR2) ) 

= a'P(Decision  at  stage  1).  (2.5) 

Since  this  Bonferroni  bound  may  also  be  overly  conservative, 
and  since  a'  will  be  less  than  a,  it  is  apparent  that  A 
will  be  small  if  it  is  positive.  If  A is  negative  the 
resulting  confidence  interval  would  be  less  conservative 
than  under  the  assumption  A = 0.  Thus,  if  we  assume 
A = 0 and  then  bound  p*  under  this  assumption,  the 
resulting  confidence  interval  will  be  such  that  P(0el)  _> 

1 - a - A,  where  A is  bounded  above  by  ( 2.5  ) . We 
proceed  with  the  assumption  that  A = 0. 
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The  expression  ( 2.4  ) for  p*  now  becomes 

p*  = P(YeR1,YeR2)  - P(YeR  ) (l-o' ) . (2.6) 

To  establish  an  upper  bound  on  (2.6  ) we  again  consider 
the  relationship  of  R^  with  R^  as  9 changes.  Our 
observations  will  pivot  on  the  relationship  of  R^  with  R2 
when  9 is  such  that 


9 


0V9 


i.e. , when  the  left  endpoints  of  R^  and  R2  coincide.  Here 
either  R2  c R^  or  R^£  R2  • We  will  designate  the  situation 
R2~  r^  as  Case  I,  and  Case  II  will  correspond  to  R^E  R2* 
Case  I and  Case  II  are  depicted  in  Figure  2.3  . In  a 
specific  application,  whether  Case  I or  Case  II  applies 
will  depend  on  the  test  parameters  involved  and  the  value 
of  a'.  Each  case  is  considered  separately. 

Case  I:  R2—  Ri  when  their  left  endpoints  coincide. 
We  examine  the  nature  of  the  relationship  of  R^  with  R2  for 
different  9 and  its  effect  on  the  problem  of  placing  a 
bound  on  p*. 

I.  a)  9 is  such  that  R^  fl  R2  = 0* 

Referring  to  ( 2.6  ) it  is  clear  that  p*  _<  0 
for  these  values  of  9. 
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Figure  2.3.  Relationship  of  R..  with  R when  0 is 
such  that  their  left  endpoints  coincide.  2 

I.b)  0 is  such  that  R1  fl  R^  f 0 and  ^ H R2  ? 0. 
Let  R1  fl  R2  = R5  and  ^ fl  R2  s r4  . with  this 
notation  we  have 

p*  = P(YeR1,YeR2)  - (1-a' ) P (YeR2) 

= P(YeR5)  - (l-a1)  (P(YeR4)  + P(YeR5)] 

= a'P(YeR5)  - (l-a')P(YeR4) . (2.7) 

Expression  ( 2.7  ) is  maximum  when  R^  = 0 or, 
equivalently  R2c  R^ . 
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It  is  clear  that,  to  establish  an  upper  bound  on  p*  for 

Case  I,  it  suffices  to  consider  only  those  0 for  which  R2 

£ . If  R2E.  R1  expression  ( 2.6  ) reduces  to  p*  = 

a'P(YeR2).  Note  first  that  if  R2—  Ri ' then  P(YeR)  _< 

P(YeR^)  = 1 - a'.  So,  one  bound  on  p*  under  Case  I is 

a'(l-a').  A smaller  bound  may  be  found  if  we  determine  the 

maximum  value  of  P(YeR2),  subject  to  R2  S.  R^  • 

For  that  purpose,  we  find  the  value  of  0,  say  0*, 

for  which  P(YeR2)  is  maximum  for  given  0g,d^  and  d2 . 

2 

Since  Y ^ Y-.  we  can  write 
A2r^ 


u(0) 

P(YeR2|0)  = J f ( t)  dt,  where  H(0)  = 0^/0, 
«■  ( 0) 


r.-l  -t/2 
t 1 e 

u(0)  = 0 d /0,  and  f(t)  = r 

02  2 1 r(ri) 


Then 


3P(YeR2 | 0) 
30 


f (u(0)  ) 


3u(0) 

30 


f(M0))  . 


Setting 


3P(YeR2  | 6) 


30 


0 
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and  making  the  appropriate  substitutions,  we  obtain 


r -1  -0nd 

f 0 d ^ — 0 - 

, 0 2 20 

d2  ~ 


60d1 


1 9 


V1  _0odi 


Finally 


U\  1 (dl"d2>  = 1. 


e*  = 


90(d2-dl) 

2r1(ln  d2  - In  d.^ 


So,  P(YeR2|0)  increases  to  a maximum  at  0*  and  then 
decreases.  Now  let 


R2* 


pdl  e0d2 

0*  ' 0* 


If  R2*  — R1'  an  uPPer  bound  on  P*  is  Clearly  a’P(YeR2*) 


under  Case  I. 


Since 


2 2 
R.  X n • r Xa I 

1 l-y-,2^ 
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and  a'  is  as  yet  unknown,  it  appears  that  it  is  not  possible 
to  determine  if  . We  note  however  that  a'  < a and 

consequently 


V 


c 


implies  that 


R2*  — 


Thus  in  most  cases,  using  the  desired  1 - a coefficient, 
it  can  be  determined  if  R2*C  R^ . When  this  relationship 
exists  we  find  that  P(0el)  > 1 - a’  - a'P(YeR2*).  So,  if 
a'  = a/ [1+P  (YeR^*)  j , we  have  P(0el)  _>  1 - a.  We  also  note 
that 


a 

1 + P(Y£R2*) 


> 


a 

2 


insures  that  the  interval  with  a'  = a/ [1+P (YeR2*) ] will 
be  less  conservative  than  the  original  intervals  in 
Theorem  2.1  using  a'  = a/2. 

"Computations  using  actual  test  parameters  given 
by  Bulgren  and  Hewett  (1973)  suggest  that  R2*  is  quite 
commonly  contained  in  R^  for  smaller  values  of  a,  say 
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ot  <_  .2.  The  improvement  over  the  intervals  of  Theorem 
2.1  will  be  discussed  later  in  this  section. 

If  ^2*5^  Ri ' t*1G  Pr°blem  becomes  more  complex. 

Here  PfYeR^fS),  subject  to  the  constraint  R2  ?=  Ri ' reaches 
its  maximum  value  when  the  left  endpoints  of  R^  and  R2 
coincide,  or  when  the  right  endpoints  coincide,  depending 
on  the  location  of  R2*  with  respect  to  R^ . Specifically, 
when  R^  3 R^ , 


^ 2 

d^  l-§-,  2rx 


® n d-i  2 

f(t)dt,if  < x2 

1 7 r 

1 2 ,£rl 


max  P (YeR. | 0) =< 


l.  a ’ _ 

1_~'2rl 


2r 

2 ,2rl 

r 


® n^i  2 

f(t)  dt  , if  > X 


i— , 2rx 


ZL  x2 

d2  “1  2r 

t 2 , ^r^ 


We  need  to  solve  the  equation 

a = a'  + a' (max  P(YeR2|6)) 
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for  a'.  The  difficulty  involved  is  the  manner  in  which  a' 
appears  in  the  limits  of  integration  in  max  P(YeR2|9). 
Consequently,  if  it  is  found  that  R2*C  R^  cannot  be 
established,  the  use  of  a'(l-a')  as  an  upper  bound  on  p*  is 
advised.  We  next  treat  the  corresponding  problem  under 
Case  II. 

Case  II:  When  9 is  such  that  the  left  endpoints 
of  and  coincide,  R2*  arguments  essentially 

the  same  as  those  used  for  Case  I,  it  suffices  to  consider 
only  the  set  of  0 for  which  R-^El  ^2'  fn  finding  an  upper 
bound  on  p* . 

If  Rx  c r2, 

p*  = P(YeR2,YeR1)  - (1-a ' ) P (YeRj) 

= PtYeR^  - (1-a')  P(YeR2) 

= (1-a*)  - (l-a')P(YeR2) 

= (1-a1) [l-P(YeR2) ] . 

Now,  since  P(YtR2)  P(YtR^)  = 1 - a',  we  again  find,  as 

we  did  in  Case  I,  that  p*  _<  a' (1-a').  This  bound  is 

achieved  by  p*  if  R^  = R2> 

We  have  shown  that,  under  the  assumption  A = 0, 

an  upper  bound  on  p*  is  a' (1-a'),  and  P(0el)  _>  1 - a'  " 

a' (1-a').  If  we  want  P(0el)  _>  1 - a,  then  set  a'  = 1 - 
1/2 

(1-a)  . In  Theorem  2.1  we  used  a'  = a/2  which  is 
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■* 

4t 


1/2 

greater  than  1 - (1-a)  . However,  for  small  a,  a/2  ^ 1 - 

1/2 

(1-a)  . Thus,  the  improvement  over  the  results  in 

Theorem  2.1  is  slight.  A further  improvement  .s 
possible  when  R^* c . Then  we  have  a*  = a/ [1+P (YeR^*) ] . 
The  degree  of  improvement  in  this  situation  over  the 
intervals  of  Theorem  2.1  depends  on  the  value  of 
P(YeR2*) . 

To  get  some  idea  of  how  conservative  the  interval 
in  Theorem  2.1  is,  some  computer  studies  were  carried  out 
using  specific  two-stage  test  parameters  given  by  Bulgren 
and  Hewett  (1973).  For  each  test  used,  the  interval 


0„d„ 
0 2 


(see  Figure  2.2  ) was  divided  by  twenty  points.  Using 
a'  = P(0 el)  was  computed  at  each  of  the  twenty  points. 


Numerical  integration  techniques  were  used  to  find 
P (YeR2 , ZeR^) . Table  2.1  gives  the  minimum  and 
maximum  values  of  P(0el)  over  the  twenty  points.  Many 
other  tests  were  included  in  this  study  but,  since 
results  on  those  were  nearly  the  same  as  for  the  tests 
given  in  Table  2.1  , their  results  have  not  been  included 

in  the  Table. 
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Table  2.2  reflects  the  improvement  which  is 

possible  if  it  happens  that  R2*E:  and  the  value  used  for 

a'  is  a/l+P  (YtR2*)  rather  than  l-(l-a)1//2.  The  desired 

confidence  coefficent  in  Table  2.2  is  .90.  As  in 

Table  2 . 1 the  minimum  and  maximum  observed  values  of 

1 /2 

P(0el)  are  given,  first  using  a'  = l-(l-a)  , and  then 

using  a'  = a/l+P  ( Ye.  R^* ) . The  confidence  intervals  using 
the  latter  expression  for  a'  are  less  conservative,  but 
the  improvement  is  hardly  dramatic.  We  do  not  sacrifice 
much  by  using  the  simpler  intervals  of  Theorem  2.1. 


Ill . A One-Sided  Confidence  Interval 

If  a one-sided  confidence  interval  for  0 is 
desired  following  the  two-stage  test,  the  same  approach  as 
used  in  Section  II  will  yield  analogous  results.  The 
confidence  interval  is  found  according  to  the  following 
rule . 

When  the  test  decision  is  made  at  stage  one,  use 


2T 


Xl\2rx 


= I.,  say. 


When  the  test  decision  comes  at  stage  two,  use 


(3.1) 
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Xo' ,2r- 


I2*  say, 
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TABLE  2.1 

Extreme  Values  of  P(9cl)  Using  Theorem  2.1 


Test  Parameters  Nominal  Confidence  Coefficients 


V61 

a 6 

.99 

.95 

.90 

. 80 

.9973 

.9855 

.9692 

.9328 

2 

.05  .05 

.9931 

.9650 

.9285 

.8497 

.9973 

.9856 

.9693 

.9  328 

3 

.05  .05 

.9931 

.9648 

.9278 

.8466 

.9974 

.9863 

.9710 

. 9365 

4 

.05  .05 

.9926 

.9600 

.9147 

. 8327 

Note  : 

The  first  entry 

in  each  cell 

is  the 

maximum 

observed 

P(0el)  over  the 

twenty  points 

used . 

The  second  entry  is  the  minimum  observed  P(8el). 


MRMi 


ut|(N  ro  (N  iT)|CN 
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TABLE  2.2 

Extreme  Values  of  P(0el)  Using  Two  Methods 


Test  Parameters 


Minimum 

P(0cl) 


Maximum 

P(0el) 


\ 


.05  .05 


.05  .05 


.10  .10 


.10  .10 


3 


.10 


.10 


.0513 

.0592 

.0513 

.0595 

.0513 

.0608 

.0513 

.0583 

.0513 

.0566 


.9272 

.9143 

.9278 

.9152 

.9290 

.9166 

.9257 

.9125 

.9241 

.9105 


.9698 

.9642 

.9636 

.9690 

.9690 

.9634 

.9701 

.9646 

.9702 

.9647 


Note:  The  first  entry  in  a cell  corresponds  to 

o’  = 1 - (l-a)l/2. 

The  second  entry  corresponds  to  a'  = a/ [1+P (YeR2*) ] . 


Y*'  •->  ♦f*- 

I ' 
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Again,  P ( 0 e I ) = 1 - a'  - p* , where 


p*  = P (D, 0el x ] - P [5, 0el 2] 


We  use  the  transformations 


and  also  define 


2T  2T 

rl  r3 

y = and  Z = _ , 


R1  “ t0/  V ,2r^  ' 


Vi  eod2 

R2  - _ 0 ' 0 _ ' 


R3  = (0,  Xa. >2r31 * 


Now  we  may  write 


p*  = P(YeR2,YeR1)  - P(YeR2,ZeR3) 


As  with  the  two-sided  intervals,  we  can  make  some 
general  observations  about  P(0el). 
i)  If 


then 


e < 9°dl 

0 < — x , 


cct'  ,2r. 


60dl 


2 

xct'  , 2r 3 


> 


2 

xa' ,2rx* 


Here  P (YeR2 , YeR^)  = P(YeR2,ZeR3)  = 0,  and  P(0el) 
ii)  If 


then 


and 


%di 


V , 2r . 


90dl 


1 9 i -T- 

V »2r. 


0od] 


- xa'/2r1 


ri  n r2  = 0. 

Since  P (YeRx ,YeR2)  = 0,  P(0el)  > 1 - a'. 
iii)  If 


then 


0Od2 

0 — 2 r 

xa',2c1 


60d2  « 

0 - 


2 

Xa,f2r1 


and 
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Furthermore , 

R ( Oe  I ) = 1 
= 1 
< 1 

We  can  utilize  the  same  approaches  to  finding  an 
upper  bound  on  p*  that  were  used  in  the  previous  section 
on  two-sided  intervals.  One  bound  on  p*  is  found  using 
the  Bonferroni  inequality,  and  the  result  is  again 
p*  a' . Consequently  we  would  set  a'  = | in  the  intervals 
( 3.1  ) to  ensure  that  P(0el)  _>  1 - a for  all  0. 

The  other  approach  is  to  set 

A = P(YeR2)P(ZeR3)  - P(YeR2,ZeR3)  = 0. 

Then  p*  becomes  P(YeR2,YeR3)  - ( 1-a' ) P (YeR2) . By  arguments 
similar  to  those  used  in  Section  II  , it  can  be  shown  that 
p*  takes  its  maximum  value  when  — Ri*  R2~  Rl'  then 

p*  = a'P(Y£R2)  <_  a (1-a1).  Or,  we  may  again  define 


-a'  - P(YeR2)  + P(YeR2,ZeR3) 
-a'  - P(YeR2) (l-P(ZeR3|YeR2) ) 

- a'  . 
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where  6*  is  that  value  of  0 for  which  P(YeR2|0)  is  maximum. 
So,  we  obtain  p*  _<  a'P(YeR2*)  when  R2*c:  \ • we  arrive  at 
the  same  terminus  as  in  Section  II  . We  may  use 
a'  = 1 - (1-a)1^2,  which  is  nearly  a/2,  or  a/ [1+P (YeR2*) ] 
if  k2*£=-  ri  anc*  a ^ess  conservative  one-sided  interval  is 
desired. 
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